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Abstract. Correlation function of complex eigenvalues oi N X N random 
matrices drawn from non-Hermitean random matrix ensemble of symplectic 
symmetry is given in terms of a quaternion determinant. Spectral properties 
of Gaussian ensembles are studied in detail in the regimes of weak and strong 
non-Hermiticity. 



1. Introduction 

Statistical ensembles of generic real, complex and quaternion matrices have first 
been introduced in the pioneering work by Ginibrc (1965) who has managed to 
derive the joint probability distribution function (j.p.d.f.) of N complex eigenvalues 
{z e } = {xi + iyi} of N x N complex {(3 — 2) and quaternion (/3 = 4) non-Hermitean 
random matrices: 

AT 

P {2 \z 1: ---,z N ) =C 2 {N) Y[\z k -z e \ 2 Y[w 2 (z e ,z e ), (1) 

N 

P$\zi, ■ • • ,z N ) = Ca{N) J] \z k - z e \ 2 \z k - z e \ 2 \zt - M 2 w 2 (ze, z t ). (2) 

k<e i=i 

Cp{N) is a normalisation constant, w 2 (z, z) is a weight function (see discussion below). 
For real matrices {f3 = 1) with no further symmetries, the reader is referred to much 
later papers by Lehmann and Sommers (1991), and also by Edelman (1997). 

Although Ginibre's derivation of Eqs. |l]) and (||) holds for random matrices with 
Gaussian distributed entries, that is for 

w 2 (z,z) = Wq(z,z) = e~ M , (3) 

we will allow the weight w 2 (z, z) to be an arbitrary benign function of z and z provided 
the normalisation Cp{N) exists. It should be emphasized [§] that such an innocent (at 
first glance) extension is quite nontrivial as it raises the question about existence of an 
underlying matrix model whose eigenvalue representation would coincide with Eqs. ([j]) 
and (||). Say, if the underlying non-Hermitean matrix commutes with its adjoint (such 
matrices are often called normal matrices, see, e.g., Oas, 1997), the interpretation of 

f E-mail address: eugene@phy.cam.ac.uk / eugenc.kanzieper@weizmann.ac.il 
§ The referee is thanked for pointing this out. 
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Eqs. (g) and (g) as a non-Gaussian j.p.d.f is correct. Notice, however, that such a 
commutativity constraint is not a must though. For example, an ensemble of weakly 
non-Hermitean matrices introduced by Fyodorov, Khoruzhenko and Sommers (1997) 
is described by j.p.d.f. of the form (|l|) and (||), see Fyodorov, Khoruzhenko and 
Sommers (1998) and Hastings (2000). 

Of particular interest is the n-point correlation function which describes a 
probability density to find n complex eigenvalues around each of the points z\, ■ ■ ■ , z n 
while positions of the remaining levels are unobserved: 

R^{z l ,...,z n ) = {N m n) t J - J d 2 Z n+1 ---d 2 Z N P ( f(z 1 ,---,z N ). (4) 

The integration measure d 2 Zg is d 2 Zg = dxgdyg. Quite often, one is also interested in 
a thermodynamic limit 

#W--,z.) = ^4« (£.-.£) (•) 

which magnifies spectrum resolution on the appropriate energy scale Sn while letting 
matrix size N tend to infinity. 

At (3 = 2, the n-point correlation function for non-Hermitean matrix model has 
also been studied by Ginibre (1965). Adopting the method of orthogonal polynomials 
introduced in the context of Hermitean random matrix theory by Mehta and Gaudin 
(1960), it is a straightforward exercise to demonstrate that i?„ (z\, ■ ■ ■ , z n ) admits the 
determinant representation 

i?£ 2) (zi, • • • ,z n ) = dct 

The scalar kernel 

K^(z,z')=w(z,z)w(z',z') Pk{z)P k {z') 

fe=0 

is expressed in terms of polynomials Pk(z) orthonormal in the complex plane z — x+iy 

J d 2 Zw 2 (z,z)P k (z)P e (z) = 5 k{ 

with respect to the measure w 2 (z, z) d 2 Z. 

For instance, the density of states and the two-point correlation function equal 

R ( ?\z) =K N (z,z) 

and 

(2,) 2 

R 2 \Zi,Z2) = K N (z 1 ,z 1 )K N (z 2 ,z 2 ) - \K N (z 1 ,z 2 )\ , 

respectively. 

Non-Hermitean random matrices at (3 = 4 have also received some attention 
in both physical and mathematical literature especially following a recent burst of 
interest to spectral properties of non- Hermitean random operators (see, e.g., Efetov 
(1997)). 

M. L. Mehta (1967) considered a non-Hermitean matrix model of symplectic 
symmetry for Ginibre's weight function Wq(z,z) = exp(— zz), and established a 
quaternion determinant structure of one- and two-point correlation functions. He 
also conjectured a similar structure to hold for all n-point correlation functions. In 
Ginibre's case, these appear in the revised 1991 edition of the Mehta's book. 



K^\zk,zi) 
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Further progress has come with development of field theoretic techniques. 
Kolesnikov and Efetov (1997), driven by possible applications in quantum 
chromodynamics (Halasz, Osborn and Verbaarschot, 1997), have formulated a 
nonlinear supersymmetry cr-model for this class of random matrices, and derived an 
expression for the eigenvalue density in a somewhat richer model [see Eq. (|2l]) below] . 

More recently, yet another field theory approach (aka replica method) was 
outlined by Nishigaki and Kamenev (2002). There, the well known Mehta's expressions 
were reproduced for one-point correlation function in case of the very same Ginibre's 
weight function Wq(z, z). Unfortunately, both mentioned techniques run into obstacles 
when one attempts to study higher-order correlation functions whilst replica cr-models 
seem to reliably provide asymptotic expansions only. 

A different route has been chosen by Hastings (2000) who suggested that there 
exists a mapping of non-Hermitean random matrices of symplectic symmetry onto a 
fermion field theory. Even though the method might have been potentially applicable 
to a study of n-point correlation functions in the bulk of a complex spectrum, these 
have not explicitly been worked out beyond the two-point correlation function. 

Our paper reports on a comprehensive treatment of integrable structure of non- 
Hermitean random matrix models at (3 — 4. It sets a transparent and coherent 
framework to study all n-point correlation functions: while easily applied to reproduce 
the results of aforementioned studies and extend them to higher-order correlation 
functions in the spectrum bulk, it may go much farther and serve as a proper starting 
point to explore eigenvalue correlations near the spectrum edges and/or address the 
issue of universality (for a recent review of the universality phenomenon in the context 
of Hermitean random matrix models see, e.g., Kanzieper and Freilikher, 1999). 

The paper is organised as follows. Section 2 announces a most general form of 
n-point correlation function whatever the weight function in Eq. (^) is. A proof 
is given in Section 3. A concept of skew orthogonal polynomials which are central 
to performing explicit calculations is elaborated in detail in Section 4. There, exact 
expressions for skew orthogonal polynomials are given in terms of multi-fold integrals. 
For the Gaussian weight Eq. (f2l|), the polynomials are evaluated explicitly. In Section 
5, n-point correlation functions for the Gaussian weight [Eq. (pl|)] are derived for 
finite N as well as in the large- N limit. Section 6 contains concluding remarks and 
briefly mentions further possible applications of the formalism developed. 

2. Correlation function at /3 = 4 and eigenvalue depletion along real axis 

For symplectic ensemble, the following representation holds for n-point correlation 
function: 



Here Q det stands for a quaternion determinant (Dyson, 1972). The 2x2 matrix 
kernel 



. . . , z n ) = Qdet K ( N '(z k ,Zi) 



J k,t=l,—,n 



(6) 




(7) 



where the 'prekernel' kn is 



2N-1 




(8) 



k,l=0 
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and M~ l is an inverse to the real antisymmetric matrix M with the entries 

M k i = / d 2 Z (z - z) w 2 (z, z) [p k (z)pi(z) - pi(z)pk(z)} . (9) 



The polynomials Pk{z) are arbitrary provided the inverse M _1 exists. Since the 
matrix M is antisymmetric, the formulas would be simplest had it contained N copies 
of the 2x2 matrix 

1 

-1 

along the main diagonal. This is achieved by letting Pk{z) be skew-orthogonal 
polynomials qk(z) in the complex domain: 

(<72fc+i,92<?)s = - <92£,92fc+i)s = r k 6iu, (10) 
(<72fc+l,<72f+l>S = (<l2k,q2t)s = 0. (11) 
The skew product (/, g) 5 is defined as Q 

</, 9)s = I d 2 Z (z - z) w 2 (z, z) [f{z)g(z) ~ f{z)g(z)\ . (12) 



With this choice in mind, the prekernel kn further simplifies to 

/ a y-^ 1 q2k+i(z) q 2 k{z') - q2k+i{z') q2k(z) . . 

k n (z,z)= > . (13) 

Tk 

k=0 

In particular, the density of states and the two-point correlation functions are 
expressed as 

R[ i] (z) = (z- z)w 2 (z,z)k n {z,z) (14) 

and 

R 2 4 \zi, z 2 ) = (zi - zi)(z 2 - z 2 ) w 2 (zi, z\) w 2 (z 2 , z 2 ) 

x k n (zi,z 1 )k n (z 2 ,z 2 ) - \k n (z!,z 2 )\ 2 + \kn(zi,z 2 )\ 2 , (15) 

respectively. 

Notice that, in accordance with our solution [Eqs. (||) and (Q)], the n-point 
correlation function universally vanishes along the real axes 3m zi = 0, 



Rn\ z !,- ■ ■ ,z n ) (x Y[ [3mz{] a , a = 2 



,,2/ 



whatever the weight function w (z,z) is. It is this specific feature of spectral 
correlations in symplectic ensembles of non-Hermitean random matrices that has 
been revealed, for n = 1, in early numerical simulations due to Halasz, Osborn and 
Verbaarschot (1997). Qualitatively, such a depletion of complex eigenvalues along the 
real axis might have been anticipated after a brief inspection of both the j.p.d.f. [Eq. 
(||)] and the definition of n-point correlation function [Eq. ([|)]. 

The results announced [Eqs. (g) - (|ilf)] will be derived in Section 3. In Section 
4, we study properties of skew-orthogonal polynomials which constitute a natural 
basis to perform actual calculations of spectral fluctuations in (3 = 4 non-Hermitean 
random matrix ensembles. The latter are addressed in Section 5, where we consider an 

+ Notice a difference from the skew orthogonality arising in the context of Hermitean random matrices 
(Mahoux and Mehta, 1991). 
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ensemble of N x N random matrices associated with a Gaussian measure. Correlation 
functions in the regimes of strong (Ginibre, 1965) and weak (Fyodorov, Khoruzhenko 
and Sommers, 1997) non-Hermiticity are explicitly given there for finite N as well as 
in the limit of infinite matrices. 



3. Derivation 



To derive a quaternion determinant representation of n-point correlation function, we 
will follow an elegant idea of Tracy and Widom (1998). These authors have introduced 
generating functional 

JV 



G[f] 



d 2 Zi 



d Zn 



fc=i 



such that the n-point correlation function R n (zi, ■ ■ ■ , z n ) defined by Eq. (Q) can be 
viewed as the coefficient of ax • • • a n in the expansion of G[f] for a particular choice 



f(z) = Ylr=i a rS 2 (z ~ z r ). Assuming that G[f] admits the representation 



G[f] = y/det(I + K N f), 

where Kjy denotes the operator with 2x2 (self-dual) matrix kernel Kn{z, z') and / 
denotes multiplication by that function, Tracy and Widom have explicitly evaluated 
the coefficient of a\ ■ ■ ■ a n and found it to be equal to the quaternion determinant in 
the r.h.s. of Eq. (||). 

Hence, in accordance with this statement (which we will name the Tracy- Widom 
theorem) one has to seek a suitable representation for G[f] with the j.p.d.f. given by 
Eq. (0). This is easy. Due to the identity 



JV 



no** 

k<£ 



n (w 



fc,/j=i 



x t) Y[(Vk ~ ye) = det 

k>e 



x e 

y k r l 



one notices that 



m 

k<i 



Zk 



zt\ 2 \z k 



JV 

z/] 2 ~ z t) = det 

e=i 



s fc-n 

zk-l 



fe=l, 
1=1,- 



fe=l, 
1=1,- 



■,2N 
,JV 



■,2JV 
,JV 



To derive the latter, we have put X£ 



in mind, the j.p.d.f. P 



(4) 
JV 



and ye = z^ in the former. With this result 



can be cast into the form 

JV 

• • ,z N ) = C 4 {N) \\_{z£ - z e )w 2 (z e ,z e ) 



x det 



e=i 
Pk-i(ze) 



fe=l, — ,2JV 
£=!,-■■, N 



Here, we have replaced the sequence of monomials {z% } by arbitrary monic 
polynomials {pk(ze)} of degree k as this leaves the value of determinant intact. If 
the polynomials p k were not monic, the normalisation prefactor G±(N) would change. 
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This representation is fairly useful due to de Bruijn's (1955) integration formula 
da(Zi) ■ ■ ■ da{Zjq) det 



fk(zi) 
9k(ze) 



fc=l,— ,2JV 

e=i,—,N 



(27V)! Pf 



da(Z) [f k (z)gi(z) - fe(z)g k {z)} 



in which da{Z) is an integration measure and 'Pf stands for pfaffian. One derives the 
following chain: 



G 2 [f] cx det 



d 2 Z (z — z) w 2 (z 7 z) 



x I 1 + f( z )} \Pk(z) Pi(z) - pe(z) p k (z)] 



det 



M k t + I d z Z (z - z)w z {z,z) 



x f(z)[Pk(z)pe(z) - p e (z)p k (z)] 



det 



>kt 



d 2 Z (z — z)w 2 (z, z) 



x f(z) [ir k (z) pt(z) - Pe(z) TTk(z)] 



(16) 



(17) 



Matrix M arisen in Eq. © is defined by Eq. (§). To derive Eq. @, we have factored 
out M on the left. Notice that Eq. (|17| ) automatically bears a proper normalisation 
G 2 [0] = 1. The polynomials ir k are 



2JV-1 



TTfe 



(z) = (M-^uPiiz). 



e=o 



The matrix appearing under the sign of determinant in Eq. (17) can be 
represented as / + AB with 

A(£, z) = f(z)(z - z) x ' 2 w{z, z) (Mz),-Mz)) , 

B{z 1 i) = (-z-zfl 2 w{z,-z)( Pt f\ 

As soon as the transposition does not affect the value of the determinant, one observes 
the identity det(/ + AB) = det(I + BA). In this case BA is the integral operator with 
matrix kernel Kn{z\, Z2)f(z2) where 

K N (z ll z 2 ) = (zi - zi) 1/2 (z 2 - z 2 ) 1/2 w(z 1 , z 1 )w(z 2 , z 2 ) 
J2iPe(zi)n e Xz 2 ) -J2£Pe(zi)ire(z 2 ) 
J2tPe( z i)Ke(z2) -J2iPe(zi)^e{z2) 

Hence, we have proven that G 2 [f] = det(I + K^f). Since the 2x2 matrix Kn is 
self-dua|| Eq. (g|) follows by virtue of Tracy- Widom theorem. This completes our 
proof. 

| Indeed, since the quaternion kjv is represented by 2 X 2 matrix 

r , / K N {z 1 ,Z 2 ) -K N {z±,Z 2 )\ 
\K N (Zl,Z2) —KN[Zi, Z2) J 
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4. Skew-orthogonal polynomials 
4-1. General weight w 2 (z,z) 

We have seen in Section 2 that skew-orthogonal polynomials defined by Eqs. ( flO| ) 
( |l2| ) represent a natural basis in which calculations become simplest. These can 
explicitly be found for a general weight function w 2 (z, z) provided the integrals below 
make sense and, in a monic normalisation, are given by the following 2n-fold integrals: 

q 2n (z) = -j- J ■ ■ ■ I d 2 Z 1 ■ ■ ■ d 2 Z n Y[(z~ z t )(z - z t ) 

n 

x \\ \ z k - zi\ 2 \z k - z t \ 2 \zt - z f ] 2 w 2 (zi,zi), 

k<l 1=1 

q2n+i{z) = -j— / ■ • ■ / d 2 Z 1 ■ ■ ■ d 2 Z n TT (z - z t )(z - zt) 

/ n \ 

x z + ^2(z k + z k ) + 



(18) 



\ k=l / 

n 

x II \ Zk ~ z ^\ z k - 2 e\ 2 II W - zt\ 2 w 2 (z e ,z e ). 

k<£ 1=1 

Here c„ is an arbitrary constant which we will set to zero, c n — whilst 



(19) 



A n = J "... J d 2 Z 1 ---d 2 Z n l[ 



\z k - Zi\ 2 \z k - zi\ 2 J| \zt - z e \ 2 w 2 (z e ,Zi). (20) 

k<l 1=1 

For similar representations of skew-orthogonal polynomials arisen in the context of 
Hermitean random matrix theory see Eynard (2001). 

To prove that Eqs. ( |l8| ) and ([l9]) obey skew-orthogonality relations [Eqs. ( |io| ) 
- (12)], it is sufficient to show that (i) (q2 n , z m )s — and (ii) (g2n+i, z m )s = for 
integer < m < 2n — 1. 

(i) Consider 



(q2n, z m ) s oc d 2 Z (z - z)w 2 (z,z) 



d 2 Z 1 ---d 2 Z n 



J\ \z k - z e \ 2 \z k ~ z ( \ 2 Y[ W - z e \ 2 w 2 (z e ,z e ) 

<£ 1=1 

n n 

z m l[(z~ z t ){z ~ zt) - z m H (z - zt)(z - z t ) 



i=i 



i=i 



As soon as 



Y\_ \z k - z e \ 2 \z k - zi\ 2 \\_{zt - z e ) JJ (z - z e )(z - zt) 
k<i 1=1 1=1 

the dual quaternion kjv is given by 

-K N (ZI,Z2) K N (z 1 ,Z2)" 



9[k n ] = ( 

\~K N (Zl,Z2) K, N (ZI,Z2) / 

Self-duality is a consequence of the equality cr^^ffc^] = 8[Kj^](jy. 
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-271-1 Jn, 



= (let 



1 Zl 



1 z„ 
1 z„ 
Vl z 



-1 *1 \ 

-271-1 -2n ' 



we derive 



(?2n, 



z > s cx 



/ 1 Zl 



x det 



n+l 

rf 2 Zi • • • d 2 Z n d 2 Z n+ i (zt - zt) w 2 (z e , zi) 
e=i 

■ z\ n \ 

■ z? n X 



< n o 
4 n o 

6 n+l z n+l 

r:2n rrm 

6 n+l ' 



1 


^1 


1 


Zn 


1 


Zn 


1 


Zn+l 


Vi 


Zn+1 




— Z. 



n + l) is irrelevant, this reduces to 



(9271, 



z )s oc 



/„ n+l 
• • • / d 2 Zi • • • d 2 Z n d 2 Z n+1 Y[ {zi ~ zi) w 2 (z e , zi) 



x det 



/I z 1 
1 zx 



1 Zn+l 
V 1 Z n +1 



r,2n 



'n+l ^n+1 ' 

The latter integrand obviously vanishes for < m < 2n thus completing the proof of 
Eq. ©. 

(ii) Consider 



(q2n+i, z m ) s oc d 2 Z (z - z)w 2 (z,z) 



d 2 Z 1 ---d 2 Z n 



J\ \zk - zi\ 2 \z k - zi\ 2 \zi - zif w 2 (z£,z e ) 
<e i=i 

/ n \ n 

z m lz + J2(z k + z k ) j JJ (z - z t ){z - z e ) 



k=l 



t=l 



Z + ^(Zfc + Zk) (z - - Zl) 

\ k=l / 1=1 



Invoking reasoning we have used in (i), this is further reduced to 



(<72n+l, Z m ) S OC 



/„ n+l 
■ I d 2 Zi ■ ■ ■ d 2 Z n d 2 Z n+ i Y[ (ze - zi) w 2 (z e , zi) 
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m \ 

1 \ 



det 



Vl 
/I 



det 



Zl 

Z n +1 
Zn+l 
Zl 
Zl 



.,2)1 



=2n 
z l 



-n+1 

7 2n 

z n+l 



.,2n 

h 

r;2n 



z 



*n+l . 
=m / 

z n+l ' 
y m+l \ 

771+ 1 



71+1 



1 Z n+1 
\ 1 Zn+l 



2n 

77+1 

■In 

77+1 



2, 



m+1 



77+1 
^777 + 1 



l + l / 
+ 1 / 



+ 



The latter trivially vanishes for < m < In — 1. This completes our proof of Eq. (|l9|). 

One may also verify that the normalisation r n in Eqs. ([l0|) and ( |Tl] ) is related to 
A n [Eq. @] as r„ = A n+1 /A n . 



4-2. Gaussian weight 

While the representations obtained above are fairly useful to study, e.g., asymptotic 
properties of general skew-orthogonal polynomials and address the issue of universality 
of eigenvalue correlations in non-Hermitean random matrix theory at (3 = 4, there is 
no need to resort to them for a simple Gaussian weight Qj] 

N ( 



w G (z, z) = exp 



1 



that we will be interested in in what follows. 

In this case, skew-orthogonal Hermite polynomials are simple: 

, T .fe+1/2 
Q2k+1 (Z ) = I I 



q2k{z) 




1=0 v ; 



21 Z 




(21) 

(22) 
(23) 



Here Hk(z) are 'conventional' Hermite polynomials 

f + OO 



nk 

H k {z) = -= 

\/7T 



dte~ % {z + ity 



orthogonal in the complex plane z with respect to the measure w 2 (z,z)cPZ (Di 
Francesco et al, 1994): 



d 2 Z Wq(z, z) H}~ z 




tt(1 -t 2 ) 1 / 2 2 k k\ 



N 



Indeed, straightforward calculation in Eqs. ( |10| ) and ( |ll|) confirms that skew- 
orthogonality is met with 



r fe = 27r(l-T) ; V 2 (l + r )i/2( 2fc + li! 



(24) 



]\[2k+2 ■ 

% This weight may be thought of as originating from the matrix model H = Hi + ivH^ , with each of 
Hula = 1,2) being drawn from statistically independent Gaussian symplectic ensembles of Hermitean 
random matrices P[Ho-] oc exp { — [iV/(l + r 2 )] Tr(if 2 )}; the parameter v 2 = (1 — r)/(l + r) [see, 
e.g., Fyodorov, Khoruzhenko and Sommers (1997)]. 
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Yet another, integral representation for q 2 k{z) holds which is more suitable for 
our purposes. To derive it we introduce the function 

ft 



1 



(2£) 



H 2t {z) 



(25) 



and notice that it satisfies the differential equation 



(l + r)—^-2rzF k (z) 
oz 



(§) 



ft 1 



(2k)\ 



2fe+l 



(Z). 



The latter is readily verified by making use of the identity H' e (z) = 2lHi-\(z) 
Integrating out, we infer 

/ _../2 \ 



F k (z) = cxp 



1 + r 



1 + t V2 



1 



(2fc)! 



dz' exp — 



1+T 



H 2k+l {z') 



where 



K 



~[2iju> ff«(o) - (-1) -^p 



Summation over £ can be performed explicitly resulting in 



1 



Ok = 



r k+l 



1 - 



2 2k + 2 k\ 



H 2k+2 {0) 



due 



Taken together with Eqs. (p3) and (J2^), this brings us to an exact integral 
representation for the even-order skew-orthogonal polynomials: 

ft 



q2k(z) = — 



exp 



Nz z 



2(1 + r) 



r k+l 



d^e 



\2J j2k^.\2^J Q 



dw exp I — 



Nw 2 



2(1 +r) 



(26) 



5. Eigenvalue correlations in (3 = 4 Gaussian ensembles 

In this section, we apply our findings to explicitly work out n-point correlation function 
for /3 = 4 non-Hermitean random matrix ensemble associated with the Gaussian weight 
Wq(z, z) [Eq. (pi])]. By letting the parameter r tend to zero, a strongly non-Hermitean 
Ginibre's ensemble is recovered. Scaling r with matrix size TV as r = 1 — a 2 /2N where 
a ~ 0(1) one accesses a regime of weak non-Hermiticity (Fyodorov, Khoruzhenko and 
Sommers, 1997, 1997a) which is known to coincide with a zero-dimensional sector of 
a supersymmetry theory of disordered systems with a direction (Efetov, 1997). 

For other papers addressing non-Hermitean Gaussian ensembles of symplectic 
symmetry by field-theoretic (cr-model) techniques see, e.g., a supersymmetry 
treatment by Kolesnikov and Efetov (1999) and a replica approach by Nishigaki and 
Kamencv (2002). Unfotunately, both techniques run into obstacles when one attempts 
to study n-point correlation function whilst replica cr-models seem to reliably provide 
asymptotic expansions only (Verbaarschot and Zirnbauer (1985), Kanzieper (2001)). 
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In accordance with Eq. (|6|), the prekernel [Eq. Cty] is the only entity needed to 
evaluate n-point correlation function. Equations (|l3), (p2[), (23) and (§3) furnish the 
desired solution 



K N {Z,Z) = — 



1/2 



2tt \1 + t 

x H 2k +i {^z 
It holds for arbitrary finite N. 



N 

1-T 

k 



3/2 



'N-l 



E 

£=Q 



Lfc=0 
1 



fe+l/2 



1 



(2* + l) 



(2*) 



z' 



-(««-> z') 



(27) 



5.2. Limit of infinite matrices: TV — > oo 

The large- iV limit is different for weakly and strongly non-Hermitean regimes. 



5.2.1. Strong non-Hermiticity. As r — * 0, the prekernel simplifies to 



K N (Z,z') 



2tt 



jv-i 

E 

fc=0 



(2fc+l)!! ^ (2*)!! 



(28) 



We are interested in a thermodynamic limit A — > oo with a blown-up energy resolution 
z i—> z/5n where Sn = (A/271-) 1 / 2 . To this end we have to evaluate 



N 



1 ( * Z ' \ /r> — 

hm —kn | 7- , -=— = V27T 



r) 3 

N 



oo 

-EE 

fe=0 £=0 . 



(zV2^) 2fc+1 (z'^) 2 
(2fc + l)!! (2l)l\ 



(z <-» z') 



(29) 



An extra power of 8n in the denominator of the l.h.s. is brought about by a prefactor 
{z- Z y/ 2 {z' - z') 1 ' 2 inEq. (§. 



Double summation in Eq. (|29| ) can be performed explicitly. Denoting 



00 /c 



(z') 



a(z,z') = EE 

fc=0 £=0 

we observe that 

dcr 1 
ir =zo- + e zz , — 
oz oz 

This suggests that we look for a(z, z') in the form 



(2fc + l)!! (2£)U 

d(J I z 
= z a - e z 



(z z') 



<j(z,z')=e^ z2+z '^A(z,z'). 



As soon as 



we obtain 



dA 

dz 



dA 
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This results in 



>n on 



lim 73-Kjv ( 7— , t— ) = Trcxp tt(z 2 + z' 2 ) erf \y/n{z — z')j 



(30) 



The latter is sufficient to evaluate all n-point correlation functions by means of Eqs. 
(0) and (0). For instance, the scaled density of states [Eqs. @ and (14)] reads 



Pi\z) = (z-z) lim -^w 2 ( 7^,7^ ) ( t->7- 



1 



JV 



)jv on 



>n °n 



8ttY 2 exp(-47ry 2 ) [ dX exp(47rr : 
Jo 



A 2 ), 



(31) 



Y = 3mz. A particular rescaling used in Eq. 

n( 4 ) 



has been chosen in such a way 



that the scaled level density p\ ' (z) approaches unity at infinity, \Y\ 



5.2.2. Weak non-Hermiticity. The regime of weak non-Hermiticity is of a 
particular interest due to its close relation to Efetov's model of disordered systems 
with a direction. We reiterate that a degree of (weak) non-Hermiticity is governed by 
a parameter r which scales with the matrix size N as r = 1 — a 2 /2N, a ~ 0(1). 

The large- N limit of the sum Eq. (|ll|) [or Eq. tj(Z%] is dominated by contributions 
of terms with k such that k/N ~ 0(1). One therefore needs the asymptotics of skew 
orthogonal polynomials qk{z) at large indices k. 

Asymptotics for odd-order skew-orthogonal polynomials <72fc+i are those of 
Hermite polynomials 7?2fc+i [Eq. (22)]. Utilising the result from standard reference 
book by Szego (1939) 

( z 



H, 



2fc+l 



\2Vk 



2 2k+1 (-l) k k\ 

-j= sm z, 

V irk 



we conclude that 



<l2k+l{z) 



2 2k + 1 (-l) k k\ t T 



/2 



'2kN 



,1/2 



V2A^/ 



sm z\ 



(32) 



Here k/N - 0(1) and zN - 0(1). 

Asymptotics for even-order skew-orthogonal polynomials q2k can be read out from 
Eq. ©. Since 



1 



we derive 

q2k(z) 



a/TTT^ ky/TTT' 
k\ 



k > 1, 



1 



fe+i 



1 



(2fc + 2)!! 



H2k-\ 




(33) 



Applying further the asymptotic formula (Szego, 1939) 



H2k 



2Vk 



2 Ik {-l) k k\ 



cos z, 



we deduce 
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Equations (^) and ( |34| ) for skew-orthogonal Hermite polynomials at k 1 make 
it now possible to evaluate the large- iV prekernel. Substituting the two equations into 
Eq. (|l3|), and replacing the sum over k by an integral we come up with 

1 ( z z '\ 71-3/2 f 1 dX a 2 A 2 • r / 

hm w k n — , — = -— / — e~ a sin tt (z - z )A . 35 

When taking the limit N — > oo, the scale 5jy has been set to 5n = Ny/z/n. 

In accordance with Eq. (||), knowledge of the scaled prekernel Eq. ( |35| ) is self- 
sufficient to have evaluated all n-point correlation functions. For instance, the density 
of states reads [Eqs. (|]) and (|lj)] 



p f\ z ) = exp(-7r 2 F 2 /a 2 ) / ^ exp(-« 2 A 2 ) sinhpTrFA]. 

2a d J A 



(36) 



6. Conclusions 



A problem of eigenvalue correlations in symplectic ensembles of non-Hermitean 
random matrices has exactly been solved by the method of orthogonal polynomials. 
In close analogy with j3 = 4 Hermitean matrix ensembles, the n-point correlation 
function is given by a quaternion determinant [Eq. of an n X n matrix whose 
entries are quaternions with an image given by 2 x 2 matrices in the form of Eq. 
(0). To evaluate the latter it is convenient (but not obligatory) to introduce a set of 
polynomials which are skew-orthogonal in the complex plane [Eqs. (|l^) and (|TT|)1. The 
skew-orthogonality set by Eq. ([l2]) represents a natural basis in which calculational 
technology is most economic. 

In Gaussian random matrix ensembles, the eigenvalue correlations are described 



by the prekernel Eq. (27) which further simplifies down to Eqs. (|3C|) and ( p5| ) for 
strong and weak non-Hermiticity, respectively. These results apply not too close to 
the spectrum edges, which may also be studied within the current framework. 

Remarkably, at j3 = 4, all n-point spectral correlation functions exhibit a peculiar 
depletion of eigenvalues along the real axis Smzj = 0, 1 < I < n, where correlations 
vanish for both arbitrary matrix size N and a probability measure w 2 (z, z). As for 
the remaining nontrivial functional dependence, we expect it to be universal as well 
once a thermodynamic limit is taken. Equations ( |l8| ) and ( |l9| ) will obviously serve as 
a proper starting point to address the universality issue in either the spectrum bulk 
or near the complex edges of the eigenvalue support. 

Acknowledgments. I wish to thank P. J. Forrester and M. B. Hastings for a 
number of useful comments and references. 
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